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Abstract 



'q A renormalization-group scheme is developed for the 3-dimensional 0(2N)- 

symmetric Ginzburg-Landau-Wilson model, which is consistent with the use 
of a 1/N expansion as a systematic method of approximation. It is motivated 
by an application to the critical properties of superconductors, reported in a 
separate paper. Within this scheme, the infrared stable fixed point controlling 
critical behaviour appears at z = 0, where z = A -1 is the inverse of the 
quartic coupling constant, and an efficient renormalization procedure consists 
in the minimal subtraction of ultraviolet divergences at z = 0. This scheme 
is implemented at next-to-leading order, and the standard results for critical 
exponents calculated by other means are recovered. An apparently novel 
O | result of this non-perturbative method of approximation is that corrections to 

scaling (or confluent singularities) do not, as in perturbative analyses, appear 
as simple power series in the variable y = zt^ v . At least in three dimensions, 
the power series are modified by powers of lny. 
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I. INTRODUCTION 

Renormalizat ion-group techniques have long been established as the standard formal 
framework for understanding critical phenomena !]-§[]• Those based on effective Hamiltoni- 
ans of the Ginzburg-Landau- Wilson type require some kind of supplementary approximation 
if concrete information is to be extracted from them, and this is normally provided by per- 
turbation theory, in the form either of the e expansion introduced originally by Wilson and 
Fisher [Q or of a fixed-dimension calculation 0. For O(N) symmetric models, a systematic 
expansion in powers of 1/N provides an alternative method of approximation ||Z| — pTOf] . While 
the methods that have been used to extract estimates of critical exponents as power series 
in 1/N are informed by renormalization-group ideas, there seems to exist no method for im- 
plementing this approximation scheme within a systematic renormalization-group analysis. 
The purpose of this paper is to describe a renormalization scheme which does this. 

As is well known, the few terms of the 1/N expansions for critical exponents that it 
has been possible to calculate explicitly do not yield accurate results for the small values 
of N that characterize physically interesting systems, so the development of a new formal 
technique is not especially useful from this point of view. Our motivation for developing this 
scheme arises from our efforts to estimate critical-point scaling functions for thermodynamic 
properties of superconductors in applied magnetic fields, in particular the specific heat, for 
which experimental data have been available for some time (see, for example [|ll]]). For this 
problem, it turns out that perturbation theory does not yield a well-controlled sequence of 
approximations. As reported in an accompanying paper [|ll|, the 1/N expansion is advan- 



tageous in this respect, and yields qualitatively encouraging results, although quantitative 
agreement with existing data is still hard to achieve. 

In outline, the renormalization-group strategy we propose works as follows. Consider, for 
example, the order-parameter susceptibility xi^Qi A , A), where t is proportional to T — T c , 
Ao is a coupling constant which can be taken as temperature-independent in the critical 
region and A is a momentum cutoff of the order of the inverse of a typical interatomic 
distance. In the limit t — > 0, we expect this susceptibility to diverge as x ~ ^ 7 > where 7 is a 
universal exponent which we would like to determine. In general, a renormalization scheme 
introduces a renormalized susceptibility xr> together with renormalized parameters t and 
A, which depend on an arbitrary length scale, say yU _1 . The fact that physics is independent 
of /j, leads to a relation of the form 

X n(t, A, A/fj) = P x (£) X R(t(£),\(£),A/(£n)) (1) 

where £ is an arbitrary number. The prefactor P x (£) and the running variables t(£) and X(£) 
depend on how the renormalization has been implemented. Now, the function xr(A A, A//x) 
diverges when its first argument, t, approaches zero, but the nature of this singularity is 
hard to determine directly. Suppose, however, that the renormalization scheme has been 
constructed in such a way that 

p x (£) ~ r^ /v * t{£) ~ tr 1/v * \{£) -* x* (2) 

in the limit £ — > 0, where 7*, v* and A* (the 'fixed-point' coupling) are constants. By setting 
£ = t u , we find in the limit t — > 



XR (t, A, A//i) ~ r^xn(h A*, A/(r»). (3) 

The exponent 7* will be equal to the true critical exponent 7 if xr(1, A*, 00) has a finite, 
nonzero value and the renormalization scheme must be designed to ensure that this is so. 
The exponent v* will then be equal to the critical exponent u, which characterizes the 
divergence of the correlation length. 

These calculations cannot be carried out exactly, so a renormalization scheme must 
be supplemented by some systematic means of approximation and its details will depend 
on this method of approximation. Perturbative methods of approximation, involving an 
expansion in powers of A, rely on the possibility of expressing exponents such as 7 as power 
series in e = 4 — d, where d is the spatial dimensionality (although in practice it proves 
possible also to implement them directly in three dimensions). Near four dimensions, the 
bare susceptibility x(t ,A ,A) and other thermodynamic functions diverge as A — > 00, so 
the principal requirement of a renormalization scheme is to remove these divergences from 
renormalized functions such as xr- One then finds that A* is of order e, so the perturbative 
analysis can be implemented self-consistently. 

For the reasons outlined above, we wish to make use of an alternative approximation 
scheme consisting of an expansion in powers of 1/iV with d = 3. A systematic means 
of obtaining the 1/N expansion is reviewed in section II. It is facilitated by an integral 
transformation of the Hubbard- Stratonovich type which expresses the partition function as 
a functional integral over an auxiliary field \l/, with an effective Hamiltonian H e s(^). This 
partition function can be calculated as a power series in 1/N by the method of steepest 
descent. Within this expansion, thermodynamic quantities are most naturally expressed as 
functions not of the temperature variable t, but rather of the exact inverse susceptibility 
to = x -1 - The temperature dependence of to is given implicitly by solving a constraint 
equation of the form to = $(to,Ao,A). At lowest order, this constraint equation simply 
locates the saddle point \l/ = —i(to — to) of the effective Hamiltonian H e R ($). 

In section III, a renormalization scheme is exhibited which manifestly extracts the cor- 
rect critical exponents at leading order. In the case of the susceptibility, the renormalized 
constraint equation t = $R,(t, A, A///) = 3>r(xS , A, A///) obeys a relation equivalent to (|J), 
but with the roles of t and Xr 1 reversed. The substantive issue, addressed in section IV, is 
how this renormalization scheme can properly be extended to higher orders. The usual per- 
turbative strategy fails at this point, because in 3 dimensions the constraint function $ and 
other thermodynamic quantities remain finite as A — > 00. For this reason, we shall actually 
omit the cutoff altogether in subsequent sections. We find, however, that the fixed-point 
coupling strength A* is infinite, and that divergences appear in the limit A — > 00. For this 
reason, we find it convenient to deal with the inverse of the renormalized coupling, z = A -1 , 
so that the fixed point appears at z* = 0. In the context of the 1/N expansion, then, the 
primary requirement of a renormalization scheme is not to remove divergences that appear 
as A — > 00, but rather to remove those that appear as z — > 0. This crucial observation is the 
first result that we wish to report. The second is a concrete renormalization scheme that 
actually does remove the divergences at z — 0. As shown in section IV, this can be achieved 
economically by means of a 'minimal subtraction' scheme which subtracts powers of \nz. 

We have, of course, checked that the critical exponents obtained from our renormalization 
scheme agree with those obtained long ago by other methods. However, our purpose in 
devising this novel renormalization-group strategy is not simply to reproduce old results by 



a different method. The application of this formalism that we have mainly in view is the 
calculation of scaling functions. In particular, the specific heat of a superconductor near its 
critical point, and in the presence of a magnetic field B, is expected to assume the scaling 
form C(T,B) rj B~ a / 2u C(tB~ 1 / 2u ) and we wish to estimate the scaling function C(x), for 
which experimental data are available. When the Ginzburg-Landau model is enlarged to 
include the magnetic field, the renormalized specific heat obeys a relation similar to ([[]), 
namely 



c K (i, z, b) = p c (£)c K (i(£), z(£), Br 2 ) ^ r^c^r^i, o, r 2 B) 



(4) 



By setting i(£) = 1 and £ = B 1 ^ 2 L, we can identify the scaling function as C(x) = 
L~ a ' 2v C^(l, 0, L~ 2 ), with L determined as a function of x by solution of the constraint 
equation, which then has the form x = L X ' U &(1, 0, L~ 2 ). Details of this calculation are pre- 



sented in [12]; here, our principal concern is to explain the renormalization-group formalism 
that makes it feasible. It happens that special considerations apply to the critical behaviour 
of the specific heat |18| and the impact of these on the renormalization-group formalism is 
discussed in section V. 

While the main purpose of this paper is to expose details of the renormalization scheme 



that will be applied in [0 , we also wish to report what we believe to be a novel feature of the 
renormalization group, revealed by the non-perturbative nature of the 1/N expansion, that 
is of some theoretical interest. The leading corrections to asymptotic critical singularities 
(sometimes described as confluent singularities) appear in the form of a scaling variable 
y = (A— X*)\T— T C \ UJU . According to conventional wisdom, the order-parameter susceptibility, 
say, has a singular part which can be expressed as x smg = \T — T^' 7 X (y) , where the scaling 
function X(y) has a power series expansion in y. In perturbation theory, this is automatically 
true, but the 1/N expansion suggests otherwise. In 3 dimensions, at least, we find that X{y) 
contains, in addition to powers of y, singular terms which at next-to-leading order are of 
the form y n In y. This is shown in section IV and further discussed, along with our other 
principal findings in section VI. 



II. THE 1/JV EXPANSION 

We consider the standard Ginzburg-Landau- Wilson theory for N complex fields, defined 
by the Hamiltonian density 

2 

(5) 



N \ / N V 

H = £ |V<Mr)| 2 + 1 Mr)\ 2 + ^ ( £ I^Wl'j 



in which to is taken to be linear in temperature (to oc T — To, where To is the mean-field 
transition temperature) and the coupling strength Ao to be temperature- independent. This 
is actually an 0(2A)-symmetric model; we have assembled its 2 A real fields into iV complex 
fields merely to facilitate the application to a Ginzburg-Landau superconductor which we 
have ultimately in view. Introducing sources ji(r) and j*(r) for the fields 4>i(r) and 4>*(r), 
we arrive at the generating functional for correlation functions 



N 

Z\j i ,f i ] = l[J D&J Dtfexp 



d 3 r 



N 



-W + £(*(r)fl(r)+tf(r)0,(r; 



i=i 



(6) 



The 1/N expansion is generated by the procedure described, for example, by Brezin et. al. 
1~3| . A Hubbard-Stratonovich transformation 



cxp 



AN 



N 



i^ E 



\i=l 



r r r N N 

Af DVexp - d 3 r—^ 2 -i^Y, 



(7) 



where M is a normalization factor whose value is of no consequence, brings the integral over 
(pi and 0* to a Gaussian form. On computing this integral, we obtain 



(8) 



(9) 



Z\ji,j*} = M /z^exp -NH cS (^) + J d 3 r J d 3 r'^2 3l (r)A(r,r';^)j*(r' 

where the effective Hamiltonian is 

H eS (^) = I d 3 r^-m 2 (r) - Tr r , r , In A(r, r'; ^) 
J Ao 

and the propagator A(r, r'\ \l/) is the solution of 

-V 2 + 1 + iV(rj\ A(r, r'; *) = 5(r - r') . (10) 

The connected correlation functions G il ... in (ri, ■ ■ -r n ) = (<f>* (ri) ■ ■ -(f)i n (r n )) c are, of course, 
obtained by differentiation of Z[ji,j*}. For example, the two-point function Gij{r,r') = 
SijG^(r — r') is given by 



G (2) (ri-r 2 ) = (^(r 1 )</» 1 (r 2 )) c 



5 2 lnZ 



<Jji(ri)<Jjf(r 2 ) 



/D«rA(n,r 2 ;*)e- 



-NH cS (V) 



j=j*=0 



J DV e- NH ^W 



Similarly, we can define a four-point function 

/ D^> A(n, r 2 ; *)A(r 3 , r 4 ; #) e -^e ff (*) 



GW(n,r 2 ;r3,r 4 ) 



f D^ e- NH ««W 



(11) 



(12) 



Owing to the factor of iV that multiplies H e g (\l/), the 1/iV expansion is now generated by 
the method of steepest descent. It proves convenient to formulate this expansion in terms 
of a parameter t , defined by 

t = r( 2 )(o), (is) 

where T^(p) is the inverse of the Fourier transform of G^ 2 \r — r'). We write 



tt(r) = -i ( to -to-j}8)+ (2iV)- 1/2 ^(r) 



(14) 



where 5 is defined by imposing the constraint (ip(r)) = 0. Since 5 has contributions of all 

orders in 1/N, we set 5 = 5 + N~ 1 5i H . The position-independent value ^o — ~i(to — to) 

locates the saddle point of H e g up to corrections of order 1/N, and the propagator A(r, r'; \I/) 
is given by A(r, r'; \&) = A(r — r') + 0(iV~ 1//2 ), where 



A(r) 



d 3 p 



D tp.r 



(15) 



(2vr)3(p2 + t ~ o) - 

The two quantities NH e ff(ty) and A(r,r'; \l/) can now be expanded in powers of 1/iV. For 
NH e s(ty), discarding irrelevant terms that are independent of ip(r), the result is 



NH. 



cir 



iiV 1 / 



1/2 V^ r 



A, 







to-f{ioAo) d 3 r^{r) 



+ - / <Prd z r'i)(r)D- l {r-r')il)(r') 



-iN- 1 ' 2 

-N- 1 



1 



d 3 r dV d 3 r" A 3 (r, r', r")ip(r)ip(r')ip(r") 



6^2 

— / d 3 r dV dV d 3 r'" A 4 (r, r', r", O^WO^'WO 
.16 J 

--S f d 3 r d 3 r' d 3 r" A 3 (r, r', r")^{r)^{r') + O (n^ 12 ^ 
where the if) propagator is given by 

D-\r - r') = \ l 5(r - r') + -A(r - r')A{r' - r) 



and the vertex functions are 
A„(ri,---,r 
The function /(to, Ao) is 



A(ri - r 2 ) A(r 2 - r 3 ) ■ ■ ■ A(r n - n) . 



to - ^o - -AoA(O) - -A 5 / d 3 rAM;r) . 



(16) 

(17) 
(18) 
(19) 



Although in principle the counterterm S is determined by the constraint (i^(r)) = while the 
definition fll3|) determines to as a function of to and Ao, we see from the expansion ( |16"D that 
practical calculations invert this logic. Thus, the term proportional to N~ 1 5o is quadratic in 
ip(r), so we arrange for fllTf ) to be true by adjusting the value of 8 . On the other hand, the 
term that is linear in ip(r) provides a counterterm that is adjusted to make (ipi^)) vanish, 
yielding a constraint (or gap) equation 

t = f(i ,\o)+0(N~ 1 ) (20) 

that implicitly determines to as a function of to and Ao- 
For the propagator A(r, r'\ *£?) we have the expansion 

A(r, r'\ V) = A(r - r') - iN- 1/2 ^-= f d 3 n A 2 (r, r'; ri)^(ri) 



-AT 



+iN- 3 ' 2 



V2- 
d 3 r ± d 3 r 2 A 3 (r, r'; n, r 2 )^(r*i)^(r 2 ) 



-S Jd 3 r 1 A 2 (r,r';r l ) 
-= J d 3 r 1 d 3 r 2 d 3 r 3 A 4 (r, r'; r h r 2 , r 3 )^(ri)^(r 2 )^(r 3 ) 



2^ 



+V25 J d 3 n d 3 r 2 A 3 (r, r'; n, r 2 )^(ri) 



+ 0(iV" 



(21) 



with 



A n (r, r'\ ri, ■ ■ ■ , r„_i) = A(r - ri)A(ri - r 2 ) • • • A(r n _i 



(22) 



The two expansions given in (O) and (21) can be summarized diagrammatically by 
the Feynman rules given schematically in figure 1 for calculating the correlation functions 
Gi 1 ...i n (ri, ■ ■ -r n ). In practice, it is most convenient to interpret these in momentum space. 
Then the and ip propagators are 



A(p) = (p 2 + t 



D(p)=\x o l + U(p) 



with 



U(p) 



d 3 k 



2 7 (27r) 3 [fc 2 + t ][(fc + p) 2 + t ] 8vr|p| 



tan 



\P\ 



,2t 



rl/2 



(23) 
(24) 

(25) 



~— ] /2 

For future use, we note that 11(0) = at , with a = 1/167T, while for large momenta we 
can write 



n(p) 



b 

\p\ 



rl/2 



t, 



-n 



p 



to 

p 2 



(26) 



with b = 1/16 and 



n(r) = ^r- 1 / 2 tan- 1 (2r 1 / 2 ) = Aa (l - - r + 



57T 



(27) 



For our immediate purposes, we require explicit expressions at next-to-leading order in 
1/N for the constraint equation (ip(r)} = 0, the two-point function G^(p) and the truncated 
four-point function T^(pi) = G^ A > (Pi) / ITt=i G^ 2 \pj), which is one-particle-irreducible with 
respect to the (f) propagator A(p). These are shown diagrammatically in figures 2, 3 and 4 
respectively. The constraint equation is 



t = f(i ,Xo) + ^A + 0(N- 2 ) 



h 



Ar 



d 3 k 



2 J (2tt) 3 
where A(i , A ) is the integral 

A(i , Ao) = 



A(k) + N~ 



] -A(i ,X )-X Q S U(0)-6 



+ o(n- 



(28) 



d 3 k d 3 k' 



■A(k) 2 A(k')D(k + k'). 



(2tt) 3 (2tt) 3 
The inverse two-point function T^ 2 \p) = G^ 2 \p)~ l is given by 

T (2) (p) =P 2 + to + N- 1 fE(p; t , Ao) - 5 1 + O (iV" 2 



(29) 



(30) 



with 



S(p; to, A ) = \ I ^ A(fc + p)D(k) 



2 J (2tt) 3 

Thus, the definition (|TJ) of t allows us to identify 

5 = £(0;t ,A ) , 
and the constraint equation is given explicitly by 



(31) 



(32) 



to — to — 



An f d 3 k 



2 J (2tt) 



A(k) + N~ 



Ao 

4 



A(£o,An)-£(0;£o,Ao)[l + A n(0)] 



+ O (N' 2 ) 



(33) 



We shall need the four-point function r( 4 )(pj) evaluated at Pi = 0. For this case, the 
multiloop integrals shown in Fig. 4 can be simplified by using the results 



d 3 k' 
(2tt) 3 



A(k'?A(k + k') 



d 3 k' 
(27p 



A(fe') A(fc + k'Y + 2A(fc') 3 A(fc + k' 



--in(fc) 

"to 

2atn" 1/2 
A; 2 + 4t 



at 



7-3/2 



1 8t 



fc 2 + 4t (A: 2 + t ) 2 



together with the value of ^ given in (|32|). We then find 

T( 4 )(0) = N-'D^O) + N- 2 [B^to, A ) + S 2 (<o, A ) + S 3 (<o, A )] + 0(iV- 3 
where the remaining integrals are 

8 3 



1 r d 3 k 

B x (t ,Xo) = -ai 1/2 D(0) 2 j j^y 3 D(k) 

r d 3 k 1 

5 2 (t ,A ) = -D(0) f—DWjp-^ 



(P + 4t 2 ) 2 (P + t )(F + 4t ) 



-S3 (to ; Ac 



d 3 A; 



(27T) 



(27T) 



1 2 



P+t 



at 1/2 £>( ) 



fc 2 + Mr 



with L>(0) = Aq +at 



-1 J_ „+-l/2 



-1 



(34) 



(35) 
(36) 

(37) 
(38) 
(39) 



III. RENORMALIZATION GROUP AT LEADING ORDER 

In the limit iV — ■> 00, the one-particle-irreducible four-point function T^(pi) vanishes, 
and so do the higher multi-point functions. The theory is therefore effectively Gaussian, all 



8 



correlations being determined by the two-point function G^{p) = (p 2 + t )~ 1 - The analysis 
of critical-point behaviour now reduces essentially to determining the dependence of to on 
the temperature-like variable to from the constraint equation ([33]). This equation can be 
written as 

t - t 0c = to + 2a\ il /2 , (40) 

where t 0c is the critical value of t , corresponding to i = 0. It represents fluctuation 
corrections to the mean-field transition temperature and is given formally by the divergent 
integral 

A f d 3 k 1 . , 

** = -f/p^F- (41) 

Near the critical temperature, the leading behaviour of t is clearly i oc (t — t 0c ) 7 , with 
7 = 2. According to standard renormalizat ion-group ideas, we should be able to find a fixed- 
point value of the coupling strength, Aq, for which this power law is exact, and this value 
is clearly Aq = oo. Within a suitable renormalization scheme, we might expect this fixed 
point to correspond to a finite value of a renormalized coupling A. The usual motivations 
for renormalization, namely the removal of ultraviolet divergences or the exponentiation 
of logarithms of to — ^oc are absent here, because the constraint equation and correlation 
functions are finite, and contain no logarithms. Nevertheless, let us define a renormalized 
temperature variable t, a renormalized coupling A and a renormalized inverse susceptibility 
fby 

A = fJ.Z x {\)\ (42) 

t = t 0c + fi 2 Z t (X)t (43) 

to = t^i (44) 

where, as usual, /i is an arbitrary renormalization scale which serves to make t, t and 
A dimensionless. We specify the renormalization factors Z t and Z\ by the normalization 
conditions 

r( 2 V = o,t = i) = / i 2 , (45) 

lim AT (4) ( Pi = 0,t=l)=fi\. (46) 

In practice, our Feynman rules yield correlation functions as functions of t rather than t, and 
these conditions are more conveniently expressed as the requirements that to ~ ^o c = Z t fi 2 
and -D(O) = /iA when t = 1. We find 

Z A (A) = (l-aA)- 1 (47) 

Z,W = —y (48) 



Expressed in terms of A and t, the constraint equation (f40"D is 



t = (1 + aX)- 1 [(1 - aX)i + 2aAt 1/2 j . (49) 



Clearly, on setting A = A* 



we do find the exact power-law behaviour t = t 1 ^. More 



generally, using the scaling fields z = A 1 — a and r = (l+z/2a)t, we can express the inverse 
susceptibility t in the scaling form 



t = t~<T{x) 



(50) 



where the scaling variable is x = zr wv /2a with up = 1 and the universal scaling function is 
given by 



T(x) = (2x)" 2 [l + 2x - VTT^x] = 1 - 2x + 5x 2 + 0(x 3 ) . 



(51) 



Formally, the scaling relation fl50| ) and the values of the exponents can be deduced from the 
renormalizat ion-group equation, which we find it convenient to formulate in terms of the 
variables t and z. Typically, given a thermodynamic quantity v4.(Ao,£o), ^ s renormalized 
counterpart A R (z,i) is given by 



A(X ,to) = Z A (z)n UA A R (z,t) 



(52) 



where Da is the canonical dimension of A. The criterion for determining the renormalization 
factor Za{z) will be discussed in the next section. The fact that A(X ,to) is independent of 
the renormalization scale \x is expressed by the renormalization-group equation 



(3{z) 



d_ 

dz 



2 - n(zj) f-~ + D A + a A (z) 



A R (z,i) = 



where 



(5{z) 
2-ri(z) 

oa{z) 



_ttfdt\ 

m dMZM 



Solution of this equation by the method of characteristics yields the relation 

A R (z,i) = P A (£)A R (z(£),t(l)) , 
where £ is an arbitrary number and z(£) is the solution of 

with the initial condition z(l) = z, while t(£) and the prefactor Pa{£) are given by 



t{£) = t£~ 2 exp 



'[t'M'')) 



P A {£) = exp / — [D A + a A (z(tj) 



(53) 

(54) 
(55) 
(56) 

(57) 
(58) 

(59) 
(60) 
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In the case at hand, taking A = to — to c , we find (3{z) = z, r](z) = 0, D t = 2 and cr t (z) 
—2a/(z + 2a). The characteristic functions are 




z(e) = z£, t(e) = tr 2 } p t ^) = i[——\ (6i) 



and it is easy to verify by substitution into ( f^?D that the relation fl5"7| ) is satisfied. To obtain 
the scaling form fl5"0|), we can choose the arbitrary parameter £ as the solution of the equation 
t(£) = 1. The neighbourhood of the critical point t = then corresponds (assuming that 
r](z) < 2) to the limit £ — > 0. In this limit, we will generally expect that 

z{£) ~ zt (62) 

i{£) ~ fr^ 2 -") (63) 

W)-^", (64) 

where u> = /3'(0), 77 = 77(0) and v = [2 + cr t (0)] _1 are the usual critical exponents. (Here, of 
course, we have u> — 1, r\ — and i> = 1.) On setting £ ~ i 1 ^ 2- ^, the relation (|57|) becomes 



t(z,i) -^'^(zi^^A) , (65) 

with 7 = (2 — 77) z/, which can be inverted to express t in the scaling form (|50|). 

IV. EXTENSION TO HIGHER ORDERS 

It is well known that the exponents obtained in the previous section are modified when 
iV is finite, and can be expressed as power series in 1/N. The series for correlation functions 
contain infrared singularities, in the form of logarithms of to, which must be exponentiated 
to obtain the correct power laws, and our aim is to find a renormalization prescription that 
will effect this exponentiation. Experience with perturbation theory and the e expansion 
suggests that application of normalization conditions such as fl45|) and ((46|) should achieve 
this, but a strategy of this kind is unsatisfactory for two reasons. At a purely practical 
level, one obtains renormalization-group functions /3(z), etc. that are cumbersome (and, 
indeed, singular) functions of z. More fundamentally, we have, so far, no reason within 
the 1/N expansion to expect that a renormalization scheme of this kind really will produce 
the correct exponentiation. We should therefore consider just what kind of renormalization 
scheme is needed. 

A simple criterion becomes apparent, if we assume that the renormalization-group anal- 
ysis will produce relations of the kind exhibited in (|5T|). On the left hand side, the function 
A(z, t) has infrared singularities when its second argument, t, approaches zero. On the right 
hand side, these singularities are removed by the condition t(£) = 1, but there remains the 
danger that they might reappear through the limit z(£) — > 0. Evidently, the singularities 
will be correctly exponentiated into the prefactor Pa{£), provided that A (z(£), 1) remains 
finite and non-zero in the limit z(£) — > 0. The primary requirement of a renormalization 
scheme is to ensure that this is so. 

To see how a renormalization scheme might work, let us examine the divergences that 
might occur in the unrenormalized theory when t is fixed to a non-zero value. The Feynman 
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diagrams generated by the 1/iV expansion are topologically similar to those in a theory of two 
fields <fi and ip, with propagators A(p) = (p 2 + i ) and D(p) = [z + n(p)]~ respectively, 
and a single interaction vertex corresponding to (f) 2 ip. Here, we use the notation z$ = Xq 1 . 
With to fixed, there are no infrared divergences. In the ultraviolet regime \p\ — > oo, we have 
A(p) ~ 1/p 2 and, provided that z ^ 0, D(p) ~ constant. Consider, then, a subintegral that 
is one-particle-irreducible with respect to both propagators, that attaches to m external (f> 
legs and n external ip legs, and has L loop integrations. It is simple to show that the 
superficial degree of ultraviolet divergence of this integral is 

^uvO?o ± 0) = 4 - L - m - In . (66) 

There are in fact only two diagrams that diverge, namely those with L = 1 and (m, n) = 
(0, 1) or (m, n) = (2, 0). These divergences are subtracted by the constraint (ip) = and by 
the additive mass renormalization corresponding to toe- The remaining theory is ultraviolet- 
finite in 3 dimensions. When zq = 0, however, the limiting behaviour of the if) propagator is 
D(p) ~ \p\. In that case, we have 

in 

n m (z = 0) = 3-2n-j. (67) 

There are now additional divergences which will appear as singularities at zq = 0, and these 
must be removed by further renormalization. Unfortunately, the freedom that we have 
to rescale the original field <fi and the parameters to — ^oc and zq does not yield a set of 
independent counterterms corresponding to the divergent subintegrals. Indeed, we can offer 
no direct proof that this rescaling will suffice to remove all the divergences. An indirect 
(and admittedly heuristic) argument suggesting that all the divergences can nevertheless 
be removed is afforded by the observation that the limit z — > is equivalent (with to c ~ 
—constant x Ao) to a limit Ao — * oo, to — > — oo, in which the model studied here should 
be equivalent to the non-linear a model whose renormalizability to all orders of the 1/N 



expansion is proved in fl4|. 



Here, we adopt the following pragmatic approach. First, we introduce a wavefunction 
renormalization factor Z$(z), so that correlation functions r^ n ^ are renormalized according 
to 

r< ft >(pi;A ,*o) = Z^z)' n / 2 Y { l\p t -zX^. (68) 

(These r^ n ^ are defined as usual by Legendre transformation of the generating functional 
lnZ[jj,j*] and would be one-particle-irreducible when calculated in perturbation theory.) 
The definition (|l~3|) then implies that to is renormalized according to 

i Q = z;\z)nH (69) 

and we introduce renormalized parameters z and t defined by 

z = Z z {z)ii- X z (70) 

t = t 0c + f , 2i ^±^Z t (z)t. (71) 
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The renormalization factors Z^(z), Z z (z) and Z t (z) must be chosen so as to make the two 
renormalized correlation functions Ir (j>j = 0) and r^ 4 ^(p = 0) and the constraint equation 
for the renormalized temperature variable t finite when z — > with t fixed to some nonzero 
value. As usual, there are many ways in which this might be achieved. Here, we propose a 
'minimal subtraction' scheme which, however, may well be susceptible of further refinement 
for specific purposes. At relative order 1/N, we show in Appendix A that the singular parts 
of the relevant Feynman diagrams are of the form (In z) x (infinite power series in z). The 
minimal way of ensuring finite limits as z — > is to subtract just the leading singular terms, 
proportional to In z. Proceeding in this way, we obtain for the renormalization factors 

Z4z) = l + ~ b S 3 lnz + 0(N- 2 ) (72) 

Z t( z ) = 1 -^k S3lnZ + ° iN ~ 2) (?3) 

Z z (z) = l + ~S 3 lnz + 0(N- 2 ) (74) 

and hence for the renormalization-group functions 

P(z) = ooz + 0{N- 2 ) (75) 

v ( z ) = 7] + 0(N~ 2 ) (76) 

2 + crt(z) = - + -^- + 0(N- 2 ) , (77) 



z + 2a 



where the exponents are given by 



-= i 4s S3+o(A " ) = i -3Sv +o(Ar " ) (78) 

in agreement with standard results (bearing in mind that our iV complex fields correspond 
to 2N real fields). The characteristic functions that were given at leading order by ( HTJ ) 
become 

z(£) = z£» , t(£) = tt-V-rt , P t (£) = e 1 '" ( z£UJ + 2a ) . ( 8 i) 

\ z + 2a J 

It is at least of formal interest to obtain the scaling form of the inverse susceptibility 
(pH|) at next-to-leading order, for the following reason. At leading order, the scaling function 
(pl|) has a simple series expansion in powers of x ~ zt uv ' . In perturbative realizations of 
the renormalization group, whether formulated as a systematic expansion in e = 4 — d or 
directly in d = 3 dimensions, the corresponding corrections to the leading scaling behaviour 
automatically appear as a power series in the scaling variable (A — A*)£^. These corrections 
have been studied in great detail, for example in [|I3|-|TT|1 . In the 1/N expansion studied here, 



this is not so. After the minimal renormalization that removes terms proportional to In z 
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from the correlation functions, there remain weaker singularities of the form z n In z, which 
will give rise to corrections of the form x n lnx. Keeping only the leading term, proportional 
to z In z, the renormalized version of the constraint equation (R3|) is 



t{z,i) 



z + 2a 



where 



2ai l/2 I 1 + — + ■ • • I + zi( 1 + -^ In z + 



- 2 ^ = ±,| _ 



(82) 



b 2 / 26 7T 2 



7T" 



(83) 



and d is a number which must be evaluated numerically, though its exact value is immaterial 
for our present purposes. The ellipsis in fl82|) indicate both terms of higher order in N~ l 
and higher powers of z. To obtain the scaling form of the inverse susceptibility t, we express 
the right-hand side of fl32"|) in the form of ( p7|) using the explicit characteristic functions ( |8~T|) 
and choose £ = £ 1 /( 2 ~ ? ?) to obtain 



t 



£1/7 



C 



c 



2a ( 1 -!- - + • • • j + zt vh f 1 + -In (^ /7 ) + ■ 



^ + 2a 
On defining the scaling field r and the scaling variable y by 

/ 1 + z/2a \ z t-'- 



(84) 



1 + c'/N + 



t 



y 



2a(l + c'/iV + ■ ■ •) 



(85) 



and the scaling function T(y) through t = T 1 T{y), we find that this scaling function is the 
solution of 



T lh [l + yT" vh 1 + — lny + — In (2aT^ /7 ) + • ■ ■ 1 = 1 
For small y, it has the expansion 

nv) 



I -71/(1 + — In j/ J H 



(86) 



(87) 



The presence of the term proportional to ?/ lny seems to be a novel result, which is discussed 
further in Section VI. 



V. RENORMALIZATION OF THE SPECIFIC HEAT 

The specific heat exponent is given in three dimensions by 

16 



a 



3u 



-1 + 



7T 2 N 



+ 0(N~ 



(88) 



As discovered long ago by Abe and Hikami [18], the fact that t a = t + 0(N x ) is asso- 
ciated with a kind of degeneracy. When exponents are obtained in the 1/N expansion by 
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exponentiating logarithms of to = t — t 0c , it is essential to account correctly for a regular 
contribution proportional to t . Specifically, writing 

C(t ) = C Q i^ a + Ci* , (89) 



it is argued in |yj that the coefficients Co and C\, considered as functions of dimension d, 
have discontinuities at d = 3 that give rise to 'anomalous' logarithms of to when the whole 
expression is expanded in powers of 1/N. Within the renormalization scheme proposed here, 
there are similar 'anomalous' logarithms of z, which prevent the specific heat from being 
multiplicatively renormalized. To be concrete, we define 



C(A , to) = N- 1 J d 3 r(<r (r)<nO)> , (90) 

where <fi 2 = X^0*0j, and find that it is given to relative order 1/N by the sum of diagrams 
shown in figure 5. We will define a renormalized function C^(z,t) by 

C(X ,to) =C(X ,0) + C 1 (Xo)io + ^- 1 Z- 2 (z)C R (z,i), (91) 



where Z t (z) = (1 + 2a/z)Z t (z) is the same renormalization factor that appears in (|7T|) . 
Dimensional analysis tells us that Ci(Ao) = CiXq 3 , where C\ is a number. Were we able to 
supply it, an analysis to all orders in 1/N of the divergences of the theory as z — > would tell 
us whether an appropriate choice of this number (at each order in 1/N) is sufficient to make 
Cr(0, t) finite. This is beyond our present skill, but we have verified that the procedure 
does work (with C\ = —A/b 2 ) at next-to leading order. In perturbative renormalization 
schemes, as is well known, an additive renormalization is also required to make a quantity 
analogous to Cr finite in the limit d — *■ 4. In that case, the additive term depends on 
the renormalization scale /x and, in consequence, the renormalized specific heat obeys an 
inhomogeneous renormalization-group equation (see, for example (|HJ). Here, by contrast, 
the first two terms on the right-hand side of ( [91]) are functions only of A and t , so Cr 
obeys a homogeneous equation of the form (|53|) . 
Finally, let us write ( |9lD in the form 

Z 2 [C(A , t ) - C(A , 0)] = C t- a + C x t , (92) 

with Cq = Cr(0, 1) and C\ = CiZ^Xq 3 . It is somewhat reassuring to find that the amplitude 
ratio can be written as 

^ = y-l + 0(iV- 1 ) (93) 



in agreement with that given in [18|], though C\ does not, in general, have a well-defined 
limit as z —>■ 0. 



VI. DISCUSSION 

We have proposed a renormalization scheme within which critical exponents and scaling 
functions for an iV-vector Ginzburg-Landau model can be estimated by means of a system- 
atic expansion in 1/N. From the field-theoretic point of view, the essential function of the 
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renormalization group is to relate the values of correlation functions in the critical region, 
where infrared singularities arise from the vanishing 'mass' t ~ T — T c , to their values when 
t is of order 1, as exhibited in (|57|). To ensure that critical singularities are correctly ex- 
ponentiated, a renormalization scheme must ensure that these singularities do not reappear 
from the running of other parameters. In the context of the 1/N expansion, this potentially 
happens when the inverse quartic coupling constant z ~ A -1 vanishes, and the essential fea- 
ture of our renormalization scheme is to ensure that renormalized correlation functions have 
finite limits when z — > 0. In 3 dimensions, the minimal way of achieving this is to subtract 
leading powers of In z, and this procedure does indeed serve to recover the standard results 
for critical exponents. 

The main purpose of this paper is to explain how the scheme works, in preparation for 
a detailed investigation of the critical properties of high-temperature superconductors to be 



described in [[L2|]. However, in the course of the formal investigation reported here, we have 
encountered a general feature that seems to have been unsuspected hitherto. It has usually 
been taken for granted that the approach to a critical point can be described by expressing 
the Hamiltonian in the form Ti = Ti* + J2i9i@i, where 7i* is an infrared stable fixed-point 
Hamiltonian and the Oi are eigenoperators of the renormalization group at this fixed point 
(see, for example []],f§]). By expanding correlation functions in powers of the coefficients 
gi, one expects to obtain corrections to the asymptotic critical singularities in the form of 
powers of the scaling variables giT Al , with exponents Aj determined by the eigenvalues of 
the Oi. In treatments based on perturbative expansions in the coupling constant A, this 
expectation is realized automatically in the case of corrections of the form (A — A*)t" ji ' 
associated with departures of A from its fixed-point value A*. Within the non-perturbative 
approximation scheme afforded by the 1/N expansion, however, this does not happen. We 
see explicitly that the fixed point corresponds to z = A -1 = 0, and that correlation functions 
do not possess power-series expansions in z. Although scaling is maintained, in the sense 
that corrections appear in terms of the scaling variable y oc zt uu , scaling functions such as 
that exhibited in (^7j) for the inverse susceptibility are not expressible as power series in y. 
At relative order 1/N, the non-analyticity of the scaling functions is logarithmic, and 
this is probably true at higher orders also. These are not, however, the same as the well- 



known logarithmic corrections that occur at the upper critical dimension d = 4 |19| . |20| . |13| 
where 4 is a marginal operator. The latter arise from a degeneracy in the renormalizat ion- 
group equations which destroy the scaling property, yielding logarithms of r rather than 
of a scaling variable such as y. In fact, our renormalization scheme is restricted to d = 3; 
we do not know in detail how to formulate a similar scheme in general dimensions. Quite 
possibly, the singularities we encountered appear, like the Abe-Hikami specific heat anomaly, 
only at special, rational values of d, and it is not clear whether they need be logarithmic 
in general. It is also far from clear whether they are of more than academic interest. In 
principle, they should presumably be present in, for example, the specific heat of 4 He near 
the lambda transition (corresponding to iV = 1), where rather precise measurements of 
critical properties have been made [pip^f . The data are reasonably consistent with the 
assumed power-law correction, but may well not be precise enough to detect a logarithmic 
factor. Unfortunately, while the non-perturbative nature of the 1/N expansion is helpful in 
indicating the presence of these logarithms, its notoriously poor convergence make it hard 
to estimate the likely sizes of the coefficients that multiply them. 
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APPENDIX A: SINGULAR PARTS OF FEYNMAN INTEGRALS 



To implement our renormalization scheme at next-to-leading order in the 1/iV expansion, 
we need to know the singularities of the integrals fl29|) and (|37D - (|39|) near z = 0. Each of 
them can be expressed in terms of integrals of the form 



<K 



d 3 k 



b m~ /m 2N 



1W^ 

fc 4 U 2 ) ' 



(Al) 



or dlx/dz, where k = \k\ and K(m 2 /k 2 ) is proportional to k 4 as k — > 0, but approaches 
a finite value as k — > oo. The singularities arise from the region of integration where k is 
large. To evaluate them, we write 



<K 



& 



dk 



1 k(zk + b) 



2 /i,2\ 



{m/k)Ii{m /k 
zk + b 



-i 



*lp- +-i 



(A2) 



where 5*3 = 1/27T 2 and 'reg' denotes contributions that are regular when z — > 0, arising here 
from the integration region < k < 1. Defining 



./ 



(Of™ 



fm\ 2 ~ I m 2N 



\k 



k 2 



*£. 



and using the series expansion f^(x) = Y, n fnH x )> we have 

^ = E/n' )j ^ + reg, 



(A3) 



(A4) 



f.» 



where 

The integral Iq is 



S3 

5 3 



oo 



dk 



m 



i k(zk + b) i+1 \k 
00 <i/c 



n-£ 



771 8z 



i /c(z/c + 6) V k J 



m\' 



5 :H 



rfA' 



'i A;(zA; + 6) 

and for n > 1, the recursion relation 

^m r°° dk /l z 

b Ji k \k zk + b 



-— - Inz + reg 
o 



(A5) 
(A6) 

(A7) 



n-l 



S^m n mz 
nb b 



(A8) 



17 



is readily solved to yield 

." ™- ] l .1 )j / z \ n-j 



S 3 m r 



■E 

j=0 



(n-j) \b 



+ - 



mz 



h = — 



S3 / mz 



In z + reg 



(A9) 



With these results in hand, we can evaluate Ik as 



i K = -—Y-f {e) (--— 

K b f-£\ Jn \ mdz, 



mz 
■—) In, 



+ reg 



]_j_]_d_ 

£ ! \ mdz 



mz 



-x'-E^!--^) /"''-- 



+ reg 






-f 1 "/!*^?* 



( » Vft ^- 



b 2 



tfi 



' m 2 z 2 ' 
b 2 



q 2 U(q 2 
mz I 



9 \ / mz 
dq I V 6 



+ reg 
+ reg 



= -?Aln z [°° dq K(q 2 ) s(q-'-^ + ^q 2 n(q 2 ) ) + reg . 
J-00 \ b b J 

The factor of In 2 can be extracted because (d/dzY(z p \nz) = \nz(d e z p /dz e ) + rej 
p > 2£. Carrying out the q integral, we finally obtain 



7-sing 
1 K 



S^lnz 



K{Q 2 



1 d 



-1 



b ~^ ' y hdQ 

where Q(z) is the solution of 

bQ + Q 2 Tl(Q 2 ) =mz. 
Solving for Q as a power series in z, we obtain 



Q 2 fl(Q 2 



Q(z) 



7-sing 



m.2 



4a 32a . 2 \ ^ ( 1 \s\ 

I — -777 mz 4 — —r~ [mz ) + (J ( (mz) I 



o 2 



6 4 



5 S , 

-— - In, 
b 



8a 96a 2 2 / , 

1 — — mz H — — — (mz) + (J \Amzy 



K{Q 2 ) . 



(A10) 
when 

(All) 

(A12) 

(A13) 
(A14) 
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FIGURES 

FIG. 1. Elements of Feynman diagrams in the \/N expansion: (a) the ip propagator; (b) 
the counterterm i"(io>^o)j ( c ) the first few terms in the expansion of the propagator A(r,r';\l/); 
(d) elementary vertex functions arising from the expansion of NH e Q. Solid lines represent the 
lowest-order (p propagator A(r — r') and filled circles represent the counterterm 8. In (c) and (d), 
the dashed lines are the legs to which ip propagators can be attached; arbitrary sequences of these 
lines and the dashed legs are allowed. 



FIG. 2. Diagrammatic representation of the constraint equation at next-to-leading order. 

FIG. 3. Diagrammatic representation of the order-parameter 2-point function at 
next-to-leading order. 

FIG. 4. Diagrammatic representation of the vertex function T^ 4 ^ at next-to-leading order. 

FIG. 5. Diagrammatic representation of the specific heat at next-to-leading order. In the first 
diagram, the cross indicates that the loop contains two (f> propagators. 
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